Phase transitions during Cyclic-Inflation and Non-gaussianity 



m 
o 

(N 
X> 

IX, 
(N 



o 

u 

CL! 
6 

CO 

03 



> 

in 

O 

m 



Tirthabir BiswatE 

Department of Physics, Loyola University, 6363 St. Charles Avenue, 
Campus Box 92, New Orleans, LA 70118, USA 

Tomi KoivistcQ 

Institute of Theoretical Astrophysics, University of Oslo, N-0315 Oslo, Norway 

Anupam MazumdaJl 

Consortium of particle physics and cosmology, Lancaster University, LAI J^YB, UK 

(Dated: February 27, 2013) 

Typically cold inflation with positive vacuum energy density dilutes all matter except the quan- 
tum fluctuations which are stretched outside the Hubble patch during inflation. However the cyclic- 
inflation scenario, where the universe undergoes many many asymmetric cycles leading to an overall 
exponential growth, operates in the presence of thermal matter. Thus, while some of the nice prop- 
erties of cold inflation is preserved, in cyclic inflationary models thermal excitations are continually 
produced and there exists a unique opportunity of imprinting interesting thermal history in CMBR. 
In particular, we will see that a phase transition can lead to observably large higher order nongaus- 
sianities as the critical temperature is approached, and that these are increasingly divergent with 
the order of the correlations. 
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I. INTRODUCTION 



It is well known that our universe can innate over 
many non-singular cycles in such a way that every cy- 
cle undergoes a slight expansion in the scale factor of the 
universe. If the universe undergoes many cycles which 
results in 50-60 e-foldings of inflation, known as cyclic 
inflation (CI) Q, then it can explain the large scale ho- 
mogeneity of the universe. For the operation of the cyclic 
inflation scenario one requires a negative potential en- 
ergy, —A, along with thermal matter, both relativistic 
and non-relativistic, which are expected in any dynami- 
cal evolution. For instance, the relativistic particles can 
come from the Standard Model degrees of freedom and 
the cold dark matter sector (which would all be relativis- 
tic at the energy scale we are interested in). Of course, 
if the universe is locked in the negative potential energy 
phase, it will keep on inflating. However, as shown in |2|, 
depending upon the nature of the potential, one can exit 
primordial inflation through a classical transition from 
negative to positive potential regions, see also for 
the study and implementation of similar mechanisms. To 
our knowledge, CI model was the first to realize a phe- 
nomenologically successful cosmology with negative po- 
tential energy , and may have interesting implications 
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for String theory and supergravity which contains nega- 
tive vacua [111 ]. 

In the present paper we will argue that the origin of 
the "relevant" seed perturbations that were created dur- 
ing the cyclic phases, i.e. roughly speaking during the 
6 — 7 e-foldings of inflation that is seen in the cosmic 
microwave background (CMB) radiation, are statistical 
thermal fluctuations in nature. In particular we will look 
at the two point correlation function which is related to 
the power spectrum, Vq, and the higher point correla- 
tion functions which are related to the parameters /nl 
and 5hl characterizing the non-gaussian features in CMB. 
The best present constraints on the /nl come from the 
WMAP7 data and for the local type of nongaussianity 
imply — 10 < / N l < 74 pjj]. For the g N L, the constraints 
are considerably weaker and both CMB [13[ and large 
scale structure [H], EH restrict it only to be less than 
about 10 5 in magnitude. 

In a companion paper [n| , we provide a comprehensive 
discussion and derivation of how statistical fluctuations 
that exist in any thermal ensemble can source fluctua- 
tions in CMB, see also We will focus on how 
in the CI scenario the observations of the CMB could 
reveal to us novel thermal features, such as phase transi- 
tion, which is impossible to detect in the traditional cold 
inflationary scenarios. In particular, we will see that if 
the scale of the negative potential energy is close to the 
critical temperature of a phase transition, then we can 
hope to find observable nongaussian signatures. Inter- 
estingly, the signal increases as one looks at higher order 
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and tachyonic instabilities may develop in this proposal JaK It is 
worth stressing that the nonsingular scenario of Refs. [l|, 01 carl 
be based on ghost- free gravity [tj. See also |1C| for a different 



attempt at nonlocal reconciliation of AdS and our universe. 
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correlation functions. So for instance, while we find fm, 
to be too small to be detected by Planck, g NL as well as 
t nl can be large enough to be observed. In addition we 
expect their amplitude to peak over a certain range of 
scales, which is a further very distinguishing prediction 
of the thermal phase transitions we explore here. 

So, let us now provide a brief summary of how the 
cyclic evolution works [T|, [H, [2(| ■ Let us begin our story 
at the expanding phase of a given cycle with the universe 
containing some thermal matter along with a negative 
cosmological constant 2 . As our universe expands, all the 
matter components dilute and is eventually canceled by 
the negative cosmological constant causing the universe 
to turnaround and start contracting. In Einstein's grav- 
ity, this leads to an eventual collapse, but, for instance, 
in an extension of General Relativity with infinite higher 
derivatives |2ll425| it is possible to resolve this cosmo- 
logical singularity, and thus the contraction phase can 
nonsingularly transition into an expansion phase thereby 
completing the cycle. An interesting fact about this class 
of infinite derivative theories is that in the ultraviolet 
regime it resolves the cosmological singularity 3 , but in 
the far infrared the theory behaves exactly as that of the 
Einstein's theory of gravity. The bounce happens at a 
scale close to the Planckian time M~ x , where M p is the 
reduced Planck mass 4 . Once the universe departs from 
the Planckian time scale, i.e. as the universe expands, we 
are far away from the bounce and the dynamics of the 
universe is governed by Einstein's theory. Since, as will 
soon become apparent, all the interesting physics that 
we are going to focus on in this paper occurs near the 
turnaround which is far away from the bounce, the de- 
tails of the bounce mechanism should not matter. Thus, 
while we do not pretend that our assumption about the 
smooth resolution of the Big Bang singularity should be 
taken lightly, at least if such a resolution mechanism via 
bounce do exist in nature around the Planck scale, we 
fully expect our analysis to be independent of the new 
physics. 

Now, interactions between relativistic and non- 
relativistic species create entropy which can only increase 
monotonically thereby breaking the periodicity of the 
cyclic evolution. In fact, this was precisely what Tol- 
man pointed out in the 1930's giving rise to Tolman's 
entropy problem [33. l35j. However, this turns out to be 
a great asset for the CI scenario. As shown in 0, [20| . 
entropy tends to increase by the same factor in every 
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FIG. 1: Cyclic Inflation Scenario: The universe starts in a 
negative energy phase undergoing asymmetric cyclic evolu- 
tion. There occurs a last cycle when the universe makes a 
transition from negative to positive energy phase followed by 
usual radiation dominated expansion. 



cycle, while the time period of the cycles remain a con- 
stant since it is governed by |A|. This means that the 
universe must grow by the same factor, say by (I + k), in 
every cycle giving rise to an overall inflationary growth, 
see Fig. Q] Most of the advantages of standard slow-roll 
inflation is straight forwardly transferred to the CI sce- 
nario, including the production of nearly scale-invariant 
density fluctuations 20] . Tolman's entropy problem also 
finds an elegant resolution if the universe is closed [36j . 
For more details of the cyclic inflationary scenario and 
how it can solve the different cosmological puzzles (flat- 
ness, horizon and homogeneity) the readers are referred 
to P, H, [20| . In this paper our primary focus would be 
the computation of the various features in the primordial 
fluctuations that are present during the CI phase. 

The paper is organized as follows. In section HH we 
review the evolution of perturbations during the cyclic 
inflation phase. The section [TTTI then focuses on the ther- 
mal fluctuations in this scenario, in particular we con- 
sider phase transitions and their impact on nongaussian- 
ity. Tensor modes are found to be insignificant in this 
context. Section[lV]brieny concludes. The technical issue 
of the transition of the modes from sub- to super-horizon 
is treated in the appendix. 



2 The generalization to a slowly evolving negative potential energy 
required to obtain the observed spectral tilt and the graceful exit 
is discussed in |2(l and 0] respectively. 

3 In general the infrared divergence of quantum fluctuations that 
could otherwise dominate over the thermal fluctuations may also 
be resolved in nonsingular cosmologies [26| , 

4 For some of the other attempts to resolve the Big Bang/Big 
Crunch singularity through bounces when a high energy scale is 
reached, see [271433] , 



II. BRIEF DISCUSSION ON CREATING 
PERTURBATIONS FROM CYCLIC INFLATION 

In Refs. 0, a new way of obtaining nearly scale- 
invariant fluctuations has been proposed in the context 
of cyclic inflationary models. The underlying mechanism 
is very similar to the standard inflationary scenarios, but 
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the implementation is very different. To understand how 
the thermal fluctuations evolve it is useful to compare the 
cosmological time scale and the physical wavelength of 
fluctuations under consideration. Naively, if the physical 
wavelength is shorter than the cosmological scale (sub- 
Hubble phase), then the fluctuations are expected to be 
sourced by thermal statistical fluctuations in the fluid 5 , 
but once the wavelengths become larger than the cosmo- 
logical time scale, thermal correlations cannot be main- 
tained and the coupled metric and fluid fluctuations are 
governed by the usual differential equations governing an 
ideal fluid. This transition is rather similar to the usual 
transition from the sub- to the super-Hubble phase in in- 
flationary cosmology, after which the amplitude freezes. 




FIG. 2: The Last Exit: The curve in the red corresponds to 
the cosmological time-scale, while the green and the yellow 
curves represent the physical wavelengths of two different co- 
moving modes. The two modes are initially in a Mixed phase 
but then they make their last exits after which they evolve 
as super-Hubble modes. While the green wavelength exits in 
the third cycle, the yellow curve has to wait for another cycle 
to make it's last exit. The two modes experiences identical 
background cosmology, the yellow curve only lagging behind 
by a cycle as compared to the green curve. 

In the traditional monotonically expanding or con- 
tracting universe the cosmological time scale is given by 
the Hubble radius. According to this reasoning all modes 
start out in a "pure sub-Hubble" phase where the phys- 
ical wavelengths of fluctuations are always smaller than 
the cosmological scale, which has a lower bound given by 
the bounce scale: 



n 



Pb 



m: 



(i) 



where pb ~ determines the total energy density at the 
time of bounce. This phase lasts until the physical wave- 
length of the given comoving mode at the bounce point 



remains smaller than rfc. With the gradual increase of the 
wavelengths though, due to the asymmetric growth in the 
cycles, the modes start to undergo a much more com- 
plicated history. The modes now start out sub-Hubble 
at the turnaround point but exits the Hubble radius in 
the contracting phase. It re-enters the sub-Hubble phase 
near the bounce only to exit the Hubble radius shortly 
after the bounce. Finally, they cross back into the sub- 
Hubble domain in the expanding branch. 

This cyclic pattern of evolution continues for a while, 
but again as the universe undergoes inflationary growth 
over many cycles, the physical wavelength keeps increas- 
ing, so that the times that a given mode spends in 
the sub-Hubble phases around the turnarounds and the 
bounces become shorter and shorter. Eventually there 
comes a cycle when the sub-Hubble phases are shorter 
than the physical wavelength of the modes. From this 
cycle on, although technically the modes do undergo sub- 
Hubble evolutions, they are ineffective in establishing 
thermal correlations, and the modes effectively remains 
super-Hubble till the rest of the duration of the CI phase. 
We call this the "last exit" . Thus the CI scenario mim- 
ics the standard inflationary mechanism of freezing the 
amplitude of perturbations after the Hubble crossing. 

Approximately it is clear that the causal time scale in 
our scenario is given by the time period of the cycle 



(2) 



5 For a detailed discussion on the subject the reader is referred to 
our companion paper [16|| . 



where (3 ~ 0(1) constant, and therefore the modes must 
be exiting when their wavelengths ~ r. In other words, 
the modes exit very near the turnaround energy density, 
in the appendix, we provide a quantitative estimate of 
this time scale. A more physical way of thinking about 
the whole process is simply to realize that the cosmolog- 
ical time scale does not correspond to the Hubble radius 
near turnarounds and bounces, but is rather controlled 
by r and ti respectively, which in turn are determined by 
the typical energy densities near these regions. The fact 
that the Hubble scale does not always correspond to the 
cosmological scale has indeed been discussed before [13] ■ 
We now come to a very important point: Suppose k\ 
is a mode which is the first to make the last exit in a 
given cycle, then once the wave number increases by a 
factor (1 + k), the fluctuations can no longer undergo 
the last exit in that given cycle because they re-enter 
the thermal sub-Hubble evolution during the following 
turnaround, see Fig. [^J They therefore, must wait for the 
next cycle for their last exit. What this means is that, if 
k is a small parameter then the last exits always occur at 
approximately the same energy densities, ~ A, near the 
turnarounds. We know from our study of perturbations 
that the energy density (or equivalently the temperature 
in our case) at the exit controls the amplitude of fluctua- 
tions that we eventually see in the CMB. This therefore 
explains why cyclic inflation model gives rise to an ap- 
proximate scale- invariant spectrum, preserving one of the 
most robust predictions of inflation theory. 
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From the above discussion, a rather distinguishing fea- 
ture of the cyclic inflation model also emerges which was 
highlighted in [20(. The fact that in a given cycle only a 
certain range of modes make their last exit means that 
there should exist small wiggles in the primordial spec- 
trum of the form 



Vm = A 2 S{k) 



(3) 



where A is strictly a constant if we have a negative cos- 
mological constant but will vary slowly if we have a V(<p) 
giving rise to an overall tilt in the spectrum, while S(k) 
is a periodic function in In k: 



S(k(l + K)) = S(k). 



(4) 



The precise shape of S(k) is determined by how the tem- 
perature rises as the modes exit in a give cycle, which 
in turn depends on the thermodynamic properties of the 
fluid. While our previous analysis provided some tanta- 
lizing hints for such oscillatory features, the Planck data 
would be able to shed more light. For us, the important 
point is that the oscillatory feature can be used in con- 
junction with the nongaussian features, that we will soon 
discuss, to differentiate our model from other inflationary 
scenarios that also produce large nongaussianities. 



is below the critical temperature. A prototype partition 
function corresponding to critical behavior such as (JS|) is 
given by 



TlnZ 1 



V 



gT 4 - A 



T 

1 



a+2 



- 1 



(6) 



The two terms in the first line correspond to the pres- 
sure coming from the usual massless degrees of freedom, 
and the negative cosmological constant that is required 
in the CI model to provide the turnarounds. The term in 
the second parentheses contain the physics of the phase 
transition; incorporating some multiplicative modulating 
functions of temperature should not change the results in 
any significant way as we are only going to be interested 
in temperatures very close to the critical temperature, 
T c . The parameter g is the number of effective massless 
degrees of freedom times (7r 2 /30), and the parameter p 
controls the contribution of the thermal density compo- 
nent undergoing the phase transition; we consider the 
regime g 3> /i. The temperature at the turnover is To, 
and it provides the overall energy scale. In terms of these 
parameters, the cosmological constant is given by 



III. STATISTICAL THERMAL FLUCTUATIONS 
& CMB 

In many extensions of the Standard model, we expect 
our universe to have undergone phase transitions at high 
temperatures, these include Grand Unified theories (su- 
persymmetric and non-supersymmctric) as well as String 
theory. Could these phase transitions have left any ob- 
servable imprint in our sky? In this section we will pro- 
vide a simple example of thermal phase transition which 
can leave detectable signals in the CMB. First the phase 
transition parameterisation will be described, then we 
will apply the results of the companion paper [16[ to 
compute the resulting bi- and trispectrum. Finally grav- 
itational wave production is discussed briefly. 



A. Phase Transition 

Phase transitions are characterized by critical expo- 
nents. For our purpose the most relevant is the one as- 
sociated with the heat capacity, commonly denoted by 
a: 



lim Cy 



T -T r 



(5) 



a = — 1 corresponds to a first order phase transition, 
where the heat capacity diverges. We are going to restrict 
ourselves to higher order transitions with a > — 1, and 
consider the case when the temperature of the last exit 



A = p(T -> 0) = -p(T -> 0) = -T 4 

a+l 



l + (l + a)p 



(7) 



and as mentioned it is negative in the CI scenario. 
Noting the thermodynamic relation 



P(T) = 

one obtains 

p = gT 4 + A 



1 dQnZ) = T dp(T) 



V d/3 



dT 



(8) 



Cy 
V 



T 

1 

T 4 T. 



a+l 



l + (l + a)J ) - J 



4gT 3 + ti^- (l + a)(2 + a) (1- 



,(9) 
(10) 



One readily checks that this is consistent with the density 
vanishing at the turnaround T = Tq and reducing to the 
cosmological term at the zero temperature. We plot the 
equation of state in Fig. [3] 



B. Non-Gaussianity 

The expressions for calculating the power spectrum 
and the nongaussianity parameters from statistical fluc- 
tuations are derived in our companion paper |16I ). Here 
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FIG. 3: The equation of state of the fluid as a function of 
temperature when T c = 0.1 M p and To = O.OlA/p. We use a = 
0.1, g — 100 and p = 0.1. At high temperatures, the radiation 
dominates and w = 1/3, at low temperatures the cosmological 
constant is important and w = — 1. The divergence occurs at 
T — Tq where the density vanishes. 



we summarize the results: 



i 



NL 



1 



jj A(T) 



1 



hp{2p' + Tp") 
24T(p'f 

25p 2 [3 (p' + Tp") + T 2 p"'} 



7 2 / A 2 {T) 



243T 2 (p>) 



where 



A(T) 



9 + 3w + 2r 
6(1 + w) 



1 



(l + w)p(2p' + Tp") 



Tp' 2 

- and 7 = 2%/2tt 3/4 ps 6.7 
P 



(11) 
(12) 

(13) 

(14) 

(15) 
(16) 



We have used the value of 7 that corresponds to a Gaus- 
sian window function [l6j], and have set = 1 since © 
already incorporated the energy coming from the cosmo- 
logical constant, radiation and the thermal matter un- 
dergoing phase transition . 

Before trying to compute and plot these different quan- 
tities, let us try to understand the physics with the help 
of some approximations. To begin with we have 6 param- 
eters: T e , T c , To, p, g, a, where T e is the exit temperature. 
However, not all of these parameters are independent or 



In the CI scenario, one also requires the presence of non- 
relativistic species |2dl |. but it's always present in small amounts, 
and in fact, decays before the onset of the contraction phase. 
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FIG. 4: The turnover temperature To of the fluid as a function 
of A. For the red solid line, p — 0.1 and for black dashed line, 
p = 1/1000. 



important. First, we can choose g to be the number ap- 
propriate for Standard Model, g - (tt 2 /30)O(100). The 
value we use for numerics is g — 100. We note in passing 
that one of the virtues of the CI model is that there is no 
need for reheating and the universe could have only con- 
tained the Standard Model degrees of freedom all along. 
Also, our computation of the exit temperature in terms of 
To in the appendix was done for radiative matter sources 
and since p <C g, the result remains applicable to our 
case; as also seen in Fig. [31 the radiative matter domi- 
nates the system when T > T). Hence it is justified to 
fix 



T = 0.8T e 



(17) 



To consider a definite and relevant phase transition, we 
can use the critical exponent corresponding to 4> 4 theory. 
It has been computed/measured using various methods, 
and the results converge to a sa 0.1, see for instance [38j ] - 
This is the value we use for numerical examples unless 
otherwise specified. 

We are thus left with three parameters: Tq,p and T c . 
Instead of the latter we shall refer to 



A = (T e - T c )/T c 



(18) 



However, the amplitude of CMB fluctuations are known 
and this can be used to constrain one of the parameters, 
a convenient choice is Tq , which gives the overall scale of 
the cycles. We show in Fig. [4] solutions for the tempera- 
ture solved by matching the amplitude of the spectrum 
with the observed WMAP value. To conclude, our pre- 
dictions depend upon two parameters: the contribution 
of the matter subject to phase transition p, and the prox- 
imity to the critical temperature, A. 

The essential feature of our results is large higher order 
nongaussianity. They contain more and more derivatives 
of p(T) and therefore become more and more divergent 
as T approaches T c . Accordingly, the Jnl parameter 
remains typically small, as seen in Fig.[5J The parameter 
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FIG. 5: Nongaussianity }nl as a function of A. For the red 
solid line, p = 0.1 and for black dashed line, p = 1/1000. 



FIG. 6: Nongaussianity gML as a function of A and p. The 
parameter is negative and diverges as T — > T c . 



approaches a constant as T — > T c and another constant 
far away from the critical temperature. The next order 
nongaussianity <?at l is already much more interesting: the 
parameter diverges as T — > T c . To understand this, note 
that the leading term in this limit, from (fT3|). is due to 
the third derivative of the density that goes like p'" ~ 
A" -2 . However, one should also take into account the 
divergent terms in the denominator, which is given by 
A 2 ~ r 2 ~ p" 2 ~ A 2 ^""), as seen from Eqs. ([141150. The 
result is that the nongaussianity diverges as gNL ~ A~ Q . 

We illustrate this in the Fig. [5] where the nongaussian- 
ity is logarithmically plotted as a function of p and A. 
While the fact that we only get an interesting signal for 
A < 10~ 4 may appear to be a fine-tuning, that is not so. 
We remind the readers that to be consistent with the ob- 
served spectral tilt of the spectrum To and hence T e has 
to vary. For n s ~ 0.95, T e varies by approximately 30%, 
which means as long as T c lies within this range, we can 
observe a signal. In fact based on our results, it seems 
that rather than having a broad <7nl spectrum, we expect 
a sharp peak in k, which is another distinctive feature 
of the thermal signal. The higher order nongaussianities 
will peak at the wavelengths that cross the horizon as the 
phase transition occurs. 

The dependence on hte critical exponent very near the 
transition is depicted in Fig. [7J In brief, for positive 
a, whenever we are close enough to the critical temper- 
ature, the g^L is large. In general, the nongaussianity 
corresponding to the ./V-point function diverges near the 
critical temperature as (1 — T /T c ) a ^~ N \ so the higher 
order nongaussianities can be significantly larger than 
the lower order ones. This property stems from the di- 
verging behavior of the higher derivatives of the density 
p and can provide a smoking-gun signal for these models. 



C. Gravity Waves 

To close the section of CMB constraints, let us inves- 
tigate the gravitational wave production in the set-up. 




FIG. 7: Nongaussianity qnl as a function of the critical ex- 
ponent a when p = 1. For the upper line, A = 10~ and for 
the lower one, A = 10 -10 . The monotonic scaling behavior 
comes apparent only very near the transition. 



For extensive thermal systems, there are no off-diagonal 
stress energy components and as a result, there are no 
source terms for the gravitational waves. The gravity 
wave spectrum is therefore governed by the quantum vac- 
uum fluctuations similar to what happens in the standard 
scalar field driven inflationary paradigm assuming that 
the gravity is indeed quantum [39]. The gravity wave 
power spectrum is thus given by 



1 



P 



12tt 3 Mi 



(19) 



where p corresponds to the energy density of the thermal 
fluid at the time of the last exit. Therefore the tensor to 
scalar ratio is given by 



r t / s 



Vh 

Vc 



1 



3v / 3tt°7 



3^,2 



1 



p- 



T 2 M V 



dp 
df 



(20) 

Again the expression has to be evaluated at the Hubble 
crossing, k = aH . The results are plotted in Fig. [S] We 
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FIG. 8: The scalar-to-tensor ratio as a function of A. For the 
red solid line, /i = 0.1 and for black dashed line, /j, = 1/1000. 



see that the amplitude of the tensor modes is typically too 
small to be observed. For particular parameter combina- 
tions though, the scalar-to-tensor ratio diverges. How- 
ever, the cyclic inflation to occur sufficiently near these 
singular points to obtain appreciable amount of gravi- 
tational waves would seem to require considerable fine- 
tuning. Hence we regard negligible gravitational wave 
contribution a robust prediction of the scenario. 



IV. CONCLUSIONS 

We computed the different features of CMB in the 
cyclic inflation scenario by applying the formalism de- 
veloped for cosmological perturbations arising from sta- 
tistical thermal fluctuations in our companion paper [l6| . 

In particular we explored the implications of phase 
transitions during the inflationary expansion. We found 
that the nongaussianity corresponding to the Appoint 
function of thermal fluctuations diverge near the criti- 
cal temperature T c as (1 - T/T c ) a ^- N \ where a is the 
critical exponent. The higher order correlations will be 
drastically amplified over a given range of scales. This 
generic feature is an imprint quite unique to the physics 
of thermal phase transitions. Observationally the higher 
order correlations are much more difficult to constrain, 
and since they are not typically huge in cold inflation 
models, the attention has been almost exclusively upon 
the bi- and trispectrum. However, as we predict that 
the amplitude of the signal is increasingly divergent for 
higher order nongaussianities, our findings can motivate 
the study of previously neglected completely new signa- 
tures. 

Coupled with the firm prediction of negligible tensor 
amplitude and subtle wiggly features in the spectra [20| . 
the structure of nongaussianities can provide a robust 
test of the cyclic inflation scenario. 
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Appendix A: Temperature of the Last Exit 

While intuitively it is clear that once the wavelengths 
become larger than the period of the oscillations, the 
modes should essentially behave as zero modes and fol- 
low the super-Hubble evolution, in any nonsingular cyclic 
model there are always two points within a cycle, the 
bounce and the turnaround, where the Hubble radius is 
infinite. So, according to the naive definition, all modes 
are back in the sub-Hubble phase! Of course, the cru- 
cial point is that larger the wavelength of the modes, 
the shorter the time that they spent in this sub-Hubble 
phase. Very conservatively, these modes need to spend 
at least a minimum of At ~ A/c time to have any "sub- 
Hubble" evolution at all. For instance, if the time the 
mode spends is less than A/c, there is no way any thermal 
correlation could be established stretching a wavelength 
A. Thus, although technically, all modes are sub-Hubble 
during the turnaround and the bounce, that is of little 
relevance. Once a mode reaches a particular physical 
wavelength they make their "last exit" and start behav- 
ing as superHubble modes for the cycles to follow. Let 
us in fact calculate this wavelength/temperature of the 
last exit point. 

Let us denote the scale factor and temperature when a 
given mode, k, exits the Hubble radius in the contraction 
phase as a k and T k respectively: 



a k H(a k ) = k 



(Al) 



Next, let us calculate the time that this mode spends in 
the sub-Hubble phase. Ignoring any asymmetry, this is 
given by 



At 



dt 



ak da 



(A2) 



where we have chosen the convention that a = 1 corre- 
sponds to the turnaround in question, so that H(l) = 0. 
The smaller the k, the closer a k is to one, and therefore 
shorter is At, as expected. Now, in order for any ther- 
mal correlation to re-establish At > A = a/k. Clearly 
then there exists a value of k below which this condition 
can never be satisfied. The largest k for which thermal 
correlation can be reestablished is then given by 



Oik 

At = — 

k 



(A3) 



The above is an implicit equation for k via (|A1[) . Using 
the above equation, for any given thermodynamic fluid, 
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one can calculate the exit temperature, k, and a k - Let 
us now demonstrate this for radiation. 

For radiation in presence of a negative cosmological 
constant, we have 



k = a k 



A 



and At 



tan 



VA 
y/3M v 



a, 2 - a? 



(A4) 



()A3j) can now be solved numerically to obtain a e , the 
scale factor at which the last exit starts to occur. Ap- 
proximately we find 



a e k. 0.8 =*> T e w 1.25T ( 



o 



(A5) 



where To is the temperature at the turnaround point. 
More physically, all those modes whose physical wave- 
length at the turnaround point is larger than a certain 
wavelength, 



A > 



y/3M p 



1.04r 



(A6) 



they never undergo any thermal sub-Hubble evolution at 
all, they must have had their last exits in some previous 
cycle. 

One can carry out a similar analysis during the bounce, 
since the time scale of the bounce is much shorter than 
the turnaround, any mode which doesn't undergo any 
sub-Hubble phase near the turnaround also has no sub- 
Hubble evolution during the bounce. 
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